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ABSTRACT 

Spectroscopic observations of the 2006 outburst of the recurrent nova RS 
Ophiuchi at both infrared (IR) and X-ray wavelengths have shown that the blast 
wave has decelerated at a higher rate than predicted by the standard test-particle 
adiabatic shock-wave model. Here we show that the observed evolution of the 
nova remnant can be explained by the diffusive shock acceleration of particles 
at the blast wave and the subsequent escape of the highest energy ions from the 
shock region. Nonlinear particle acceleration can also account for the difference of 
shock velocities deduced from the IR and X-ray data. The maximum energy that 
accelerated electrons and protons can have achieved in few days after outburst 
is found to be as high as a few TeV. Using the semi-analytic model of nonlinear 
diffusive shock acceleration developed by Berezhko & Ellison, we show that the 
postshock temperature of the shocked gas measured with RXTE PCA and Swift 
XRT imply a relatively moderate acceleration efficiency characterized by a proton 
injection rate rj in j £ 10~ 4 . 

Subject headings: novae, cataclysmic variables - stars: individual (RS Ophiuchi) 
- acceleration of particles - cosmic rays - shock waves 



1. Introduction 

RS Ophiuchi is a symbiotic recurrent nova with various recorded eruptions in the last 
century (the last one in 1985), which erupted again recently, on 2006 February 12 (Hirosawa 
2006). RS Oph's binary system consists of a white dwarf with mass near the Chandrasekhar 
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limit, and a red giant (RG) companion star. High white dwarf mass and large accretion 
rate lead to a much shorter recurrence period of outbursts than in classical novae (where the 
donor is a main sequence star). In addition, the presence of the RG wind in the RS Oph 
system leads to the generation of an X-ray emitting blast wave that runs into a relatively 
dense circumstellar medium. 

The latest outburst of RS Oph has been observed at various wavelengths, e.g. in radio 
(O'Brien et al. 2006), IR (Monnier et al. 2006; Das et al. 2006; Evans et al. 2007) and X-rays 
(Sokoloski et al. 2006; Bode et al. 2006). The X-ray data have allowed to clearly identify 
the forward shock wave expanding into the RG wind and to estimate the time evolution of 
its velocity, v s , through the well-known relation for a test-particle strong shock: 



where k is the Boltzmann constant, T s is the measured postshock temperature and /im# is 
the mean particle mass. The X-ray emission has revealed that after an ejecta-dominated, free 
expansion stage (phase I) lasting ~6 days, the remnant rapidly evolved to display behavior 
characteristic of a shock experiencing significant radiative cooling (phase III; see Sokoloski 
et al. 2006; Bode et al. 2006). At day 6 after outburst, however, the shocked material was 
so hot, T s ~10 8 K, that its cooling by radiative losses was probably not important for the 
dynamics of the shock. Thus, the lack or the very short duration of an adiabatic, Sedov- 
Taylor phase (phase II) differs from the remnant evolution model developed by Bode & Kahn 



(1985), O'Brien & Kahn (1987), and O'Brien et al. (1992) after the 1985 outburst of RS 



The time-dependence of shock velocity has also been measured by IR spectroscopy, 
using the observed narrowing of strong coronal emission lines (Das et al. 2006; Evans et al. 
2007). Although the general behavior of the shock evolution was found to be consistent with 
that deduced from the X-ray emission, the shock velocities determined from the IR data are 
significantly greater than those obtained using eq. (1) together with the X-ray measurements 



In this Letter, we show that production of nonthermal particles at the forward shock 
through the first-order Fermi acceleration process can be deduced from these observational 
data. Several observations in the solar system and beyond show that diffusive acceleration of 
particles can be efficient in collisionless shocks and the backpressure from the energetic ions 
can strongly modify the shock structure (e.g. Jones & Ellison 1991). In particular, eq. (1) 
is known to underestimate shock velocities when particle acceleration is efficient, because 
the postshock temperature can be much lower than the test-particle value (Decourchelle et 
al. 2000; Ellison et al. 2007). We start in § 2 with a simple description of the dynamical 
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Oph. 



of T s (see Fig. 1). 
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evolution of the blast wave based on IR and X-ray observations, we then calculate in § 3 the 
maximum possible energy of particles accelerated in the nova remnant and determine in § 4 
the properties of the cosmic-ray modified shock including the energy carried off by particles 
escaping the shock region. Our conclusions follow in § 5. 

2. Shock wave evolution 

The FWZI (full width at zero intensity) of IR emission lines of coronal origin should 
provide a good measurement of the shock velocity (Das et al. 2006; Evans et al. 2007). 
Radio imaging (O'Brien et al. 2006) and IR interferometric observations (Monnier et al. 
2006) of the RS Oph remnant have shown departures from spherical symmetry. But given 
the intermediate angle between the symmetry axis of the observed bipolar structure and the 
line of sight, #~50-60 o (O'Brien et al. 2006), the largest blueshifted and redshifted velocities 
measured in the FWZI should be close to the mean expansion speed of the blast wave. The IR 
data can be modeled at first approximation by (dashed line in Fig. 1) v s (t) = A300r av km s -1 , 
where t is the time after outburst, r—tjt\ with ti=6 days, and 0:^=0 (-0.5) for t<t\ (t>ti). 
The decay of v s as t~ a5 for t>ti can be expected from a well-cooled shock. We see that 
the velocities deduced from the X-ray data are consistent with such a decay (dotted line in 
Fig. 1), although they are ~1.7 times lower than the velocities determined from the IR lines. 

The radius of the shock front, which for simplicity we assume to be spherical, is easily 
obtained by integration of v s (t): r s (t) = 2.23 x 10 14 [(1 - 2a v )r av+1 + 2a v \ cm. We do not 
consider the earliest phase of the outburst, when the shock wave traversed the binary system. 
Given the binary separation of ~1.5 AU (Fekel et al. 2000), the free expansion of the ejecta 
into the unperturbed RG wind started at t~t =l day. The outer radius of the RG wind is 
^out=MRG^, with mrg— 10-20 km s -1 the terminal speed of the RG wind and At=21.04 yr 
the elapsed time between the 1985 and 2006 outbursts. The outer radius was reached by the 
forward shock at t 2 = 24-72 days after the 2006 outburst. 

The density of the RG wind as a function of radius is given by pw{ r ) — ^RG/(47rr 2 iiRG)> 
where Mrg is the RG mass-loss rate. It can be estimated from the X-ray photoelectric 
absorption measured with Swift XRT (Bode et al. 2006) . We use equation (4) of Bode et 
al. (2006) to fit the measured absorbing column density. We do not take into account in this 
fit the data taken at t=3. 17 days, because the early X-ray emission could partially originate 
from the reverse shock running into the ejecta. We obtain MrgA*rg~4x10 13 g cm' 1 , which 
is ~5 times higher than the value estimated by O'Brien et al. (1992) for the 1985 outburst. 
Part of the difference is due to the larger ejecta velocity considered here. 
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Relatively large magnetic field of stellar origin is expected to pre-exist in the RG wind. 
Assuming that turbulent motions in the wind amplify the magnetic field By/ up to the 
equipartition value (Bode & Kahn 1985), we have B w = (87ip w kT w / pm H )°- b , where T w is 
the wind temperature, which we assume to be uniform throughout the wind and equal to 
10 4 K. In the vicinity of the shock front, further magnetic field amplification is expected to 
occur due to various interactions between accelerated particles and plasma waves (e.g. Lucek 
& Bell 2000). Assuming a time-independent amplification factor a# (expected to be of the 
order of a few), the magnetic field just ahead of the shock, B =o:bB w , can be evaluated 
using the above relations for the RG wind density and shock radius: 

B (t) = 0.047a B [(l - 2a v )r av+1 + 2a,,]" 1 G. (2) 

Given these estimates of v s , r s , Po=pw( r s), and B Q , we can now study the effects of particle 
acceleration at the blast wave. 



3. Acceleration rate and maximum particle energy 

The rate of momentum gain of nonthermal particles diffusing in the vicinity of a shock 
wave is given by (e.g. Lagage & Cesarsky 1983) 

/dp\ = p(u - u 2 ) / «o(p) + K 2 (p) \ - 1 

V dt ) acc 3 V Uo U2 ' 

where uq (112) is the upstream (downstream) component of flow speed normal to the shock 
in its rest frame and kq (/C2) is the upstream (downstream) spatial diffusion coefficient in 
the direction normal to the shock. Here and elsewhere, the subscript "0" ("2") implies 
quantities determined far upstreanQ (downstream) from the shock front, as in Berezhko 
& Ellison (1999). To estimate the spatial diffusion coefficient, k=Xv/3, we assume the 
scattering mean free path A of all particles of speed v to be A = 77 m f p r 9 (Baring et al. 1999; 
Ellison et al. 2000), where 7? m fp is a constant and r g =pc/ (QeB) is the particle gyroradius, c 
being the speed of light, Q the charge number, and — e the electronic charge. The minimum 
value T] m f p —1 corresponds to the Bohm limit. 

Because the wind magnetic field should be largely azimuthal, the acceleration efficiency 
could be significantly reduced as compared with that in a quasi-parallel shock (shock normal 
parallel to the field direction). A full treatment of this question is beyond the scope of this 



1 By far upstream, we mean ahead of the upstream shock precursor induced by the backpressure of 
energetic particles. 
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Letter (see, e.g., Kirk et al. 1996). For simplicity, we consider the shock to be quasi-parallel 
and allow for large values of ?7 m fp to account for the possible reduction of the acceleration 
efficiency. 

Given the expected field compression B2=r tot B , where r tot =u /u2 is the total compres- 
sion ratio of the shock, we have K2(j>)=kq(j>) / r\ ot . The rate of energy gain for ultrarelativistic 
particles (v =c) can then be written from eq. (3) as 

( d 4) =50.0^1-g-g o r ^ MeVs- 1 , (4) 

V at / acc r tot ?7 m fp V lCr km s 1 / 

where B is in units of Gauss. Acceleration can only occur for particles whose acceleration 
rate is lower than their energy loss rate. Concentrating here on nucleons, we can safely 
neglect the collision energy losses. Adiabatic losses, however, can limit the maximum energy 
-Etaax that particles can acquire during the shock lifetime. The rate of momentum loss due 
to adiabatic deceleration of the nonthermal particles in the expanding flow can be written as 
(dp/dfjad = pv s (r tot — l)/(3r s r tot ) (e.g. Volk & Biermann 1988), which gives for the energy 
loss rate of ultrarelativistic particles accelerated in the nova remnant 

(^) = 0.64 ^— ^ £ TeV [(l - 2a v )r + la^y 1 MeV s" 1 , (5) 
V at / ad r tot 

where -&rev is the particle kinetic energy in TeV. By equalling the adiabatic loss rate with 
the acceleration rate (eq. 4), one obtains an upper limit on the nonthermal particle energy: 

E < E niax = 68-^-a B r^ TeV. (6) 

?7mf P 

It is likely, however, that high-energy particles can escape the acceleration process before 
reaching the maximum energy given above. Following Baring et al. (1999), we assume the 
existence of an upstream free escape boundary (FEB) ahead of the shock, located at some 
constant fraction / esc of the shock radius: d-pEB = /esc?V The maximum energy that particles 
can acquire before reaching the FEB is obtained by equalling d-pEB to the upstream diffusion 
length l Q ~ Kq/v s . Using the parameters derived in § 2, one finds for ultrarelativistic particles: 

iC a ; = 136f csc -^-a B r a * TeV. (7) 

??mfp 

We see that for / csc <0.5 this size limitation of the acceleration region gives a more restrictive 
constraint on the maximum particle energy than that obtained from the adiabatic losses 
(eq. 6). 

Calculated maximum proton energies are shown in Fig. 2 for / csc =0.25 (Baring et al. 
1999) and a B /rj m{p =0.1. The quantity E^ x is the maximum proton energy caused by the 
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finite age of the shock, and has been obtained by time integration of (dE/dt) acc from t to 
t, adopting r tot = 6.5 (this value of r tot is close to what we obtain for t>ti, see Fig. 3c). 
We have slightly underestimated -E^fax f° r the initial nonrelativistic phase by assuming v—c, 
but the error is negligible given the very short duration of this phase. The production of 
proton energies above 1 TeV is a consequence of the high values of Bq (eq. 2) implied by the 
assumption of equipartition. In Fig. 2, the time for which E^ X =E^ X , i.e. the beginning 
of particle escape from the shock region is t csc =5.95 days. This time only depends on / esc , 
but not on a B /rj m{p which is a scale factor for both E^ x and E^.. It is remarkable that 
for / csc =0.25, t esc is very close to the observed transition time t±, which could explain the 
apparent lack of an adiabatic phase in the remnant evolution. 

4. Properties of the cosmic-ray modified shock 

Berezhko & Ellison (1999) have developed a relatively simple model of nonlinear diffusive 
shock acceleration, which allows to quantify the modification of the shock structure induced 
by the backreaction of energetic ions. Although the model strictly applies to plane-parallel, 
steady state shocks, it has been sucessfully used by Ellison et al. (2000) for evolving super- 
nova remnants. Given the upstream sonic and Alfven Mach numbers of the shock, which can 
be readily calculated from the parameters derived in § 2, and the maximum particle energies 
evaluated in § 3, both the thermodynamic properties of the shocked gas and the energy 
spectrum of the accelerated protons (other particle species can be neglected for evaluating 
the shock modification) are determined by an arbitrary injection parameter 77^, which is 
the fraction of total shocked protons in protons with momentum p>pi n j injected from the 
thermal pool into the diffusive shock acceleration process. We used the work of Blasi et al. 
(2005) to accurately relate the injection momentum p in j to r] in y 

Calculated temperatures of the postshock gas are shown in Fig. 3a. We see that the 
temperatures measured with RXTE/PCA and Swift/XRT can be well reproduced with 
?7i n j = 1.4xl0~ 4 and Alfven wave heating of the shock precursor. The latter process is thought 
to be an important ingredient of cosmic-ray acceleration (McKenzie & Volk 1982) and ap- 
pears to be required in this case as well to limit the shock compression ratio and acceleration 
efficiency For 77i n j = 1 5 , the test-particle approximation applies and the standard relation 
between v s and T s (eq. 1) overestimates the temperature. 

The solution shown in Fig. 3a is not unique. For example, an equally good description 
of the T s measurements can be obtained with r/ m f p =100 and r]i n j=1.9xl0~ 4 . However, all 
the solutions providing good fits to the data give about the same compression ratio and 
acceleration efficiency The latter is shown in Fig. 3b for the same input parameters as in 



- 7- 



panel (a). The two quantities plotted in this figure are (thin curves) the fraction of total 
energy flux, F Q , going into nonthermal particles, and (thick curves) the fraction of F 
escaping the shock system via diffusion of the highest energy particles across the FEEq Here, 
F = 0.5p v^ + [jg/i'fg — ^)]Pq v s, where 7 9 =5/3 is the specific heat ratio of the upstream 
thermal gas and Pq is the far upstream pressure. We see that for r]i n j=1.4xl0~ 4 and Alfven 
wave heating, e csc =10-20%, which means that accelerated particle escape is important for 
the dynamics of the shock. Compression ratios and postshock pressures calculated for this 
case are shown in Fig. 3c. 

The energy loss rate due to the escape of the highest energy particles can be estimated 
to be 

(f)„^ 4 -. 2 - F »- 2 - 4xl038 (<nl)^ 5 (8) 

for t > t csc = t%, where we have neglected the enthalpy term |/y g / (•y g — l)]P v s in the expression 
for F . For t = t\, (dE/dt) esc is ~200 times higher (for the assumed distance of 1.6 kpc; 
Hjellming et al. 1986) that the bolometric luminosity of the postshock hot plasma (Bode 
et al. 2006; Sokolovsky et al. 2006), which shows that energy loss via accelerated particle 
escape is much more effective to cool the shock than radiative losses. 



5. Conclusions 

We have shown that production of nonthermal particles by diffusive acceleration at the 
blast wave generated in the 2006 outburst of RS Oph can reconcile shock velocities deduced 
from X-ray data with velocities measured in broad IR lines of coronal origin, and account 
for the observed cooling of the shock starting as early as ~6 days after outburst. 

Using a semi-analytic model of nonlinear diffusive shock acceleration, we have con- 
strained the proton injection rate from the measured postshock temperature to be ?7i n j ^ 10~ 4 . 
We believe that the existing high-quality multiwavelength observations of this nova outburst 
could allow to further test and improve the diffusive acceleration theory. 

To our knowledge, the acceleration of particles to TeV energies in a recurrent nova rem- 
nant was not predicted. In a forthcoming paper, we will calculate the high-energy emission 
generated via interactions of this nonthermal population with the ambient medium. 



2 Although e csc has been set to zero for t < t csc , it is likely than some high-energy particles escape 
upstream from the shock at all times, because of the low scattering strength of high-momentum particles in 
self-generated turbulence (see Vladimirov et al. 2006). 
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Fig. 1. — Time-dependence of forward shock velocity as deduced from FWZI of IR emission 
lines (filled symbols: squares and circles (Das et al. 2006), triangles (Evans et al. 2007)) 
and X-ray measurements of the postshock temperature (open symbols: triangles (Sokoloski 
et al. 2006), squares (Bode et al. 2006)). The lines are simple fits to the data (see text). 
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Fig. 2. — Estimated maximum proton energy as a function of time after outburst. E^ x and 
£™ are the maximum energies caused by the finite shock age and size, respectively (see 
text). The solid curve shows the minimum of these two quantities. 
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Fig. 3. — (a) Calculated postshock temperature compared to the RXTE and Swift data (same 
symbols as in Fig. 1). Dashed curve: ?7 in j=10~ 5 . Solid and dotted curves: ?7i n j=1.4xl0~ 4 
with and without Alfven wave heating of the shock precursor, respectively. All calculations 
assume a#=2, r] m f p =20 and f esc =0.25. (b) Nonthermal energy fraction (thin curves) and 
escape energy fraction e csc [thick curves) for the same parameters as in panel (a). For 
77 in j=lCr 5 , e esc is <1CT 2 . (c) Subshock and total compression ratios (dashed curves, left axis) 
and normalized postshock pressures (right axis) in thermal (-Pth) and nonthermal (-Pcr) 
particles (-Ptot=-Pth+-PcR,) f° r ^inj=l-4x 10~ 4 and Alfven wave heating. 



